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Abstract 

Constraints on the spin-averaged Ap scattering length and effective range 
have been obtained from measurements of the pp — > pK + A reaction close to 
the production threshold by comparing model phase-space Dalitz plot occu- 
pations with experimental ones. The data fix well the position of the virtual 
bound state in the Ap system. Combining this with information from elastic 
Ap scattering measurements at slightly higher energies, together with the fact 
that the hyperdeuteron is not bound, leads to a new determination of the low 
energy Ap scattering parameters. 
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1 Introduction 



The existence of light hypernuclei, such as ^He, shows the low energy A-p inter- 
action to be strongly attractive, though not sufficient to bind the two-baryon hy- 
perdeuteron. The A-p interaction is of especial interest since it is influenced by 
the strange quark content of the A-hyperon. However, in contrast to the nucleon- 
nucleon case, due to the short lifetime of the A, direct measurements of low-energy 
A-p scattering are sparse and the resulting parameters rather poorly known. 

Bubble chamber measurements |JJ, |^, ffl, based on samples of a few hundred 
secondary events, have allowed determinations of the elastic cross section down to 
A laboratory momenta ~ 130 MeV/c. In the low energy region, where only S- waves 
are important, the spin-averaged total cross section is of the form 




Here q is the Ap centre-of-mass momentum and cl s m and r s M are, respectively, the S- 
wave scattering lengths and effective ranges in the Ap spin-singlet and triplet states. 
Separate values of these parameters have been claimed for the two spin states 
and these are shown in Fig. [|. However, the error bars are large, strongly and 
systematically correlated and hard to quantify, since such data should really only 
support the determination of an average scattering length a and effective range f 
[[J. Already for laboratory momenta ~ 300 MeV/c, the differential cross section is 
significantly non-isotropic, indicating the presence of P or higher waves P], E3], and 
so it is not surprising that the S'-wave parameters deduced from such experiments 
depend upon the upper momentum cut assumed. 

Values of the scattering length and effective range have also been deduced through 
the study of the Ap final state interaction (FSI) in the K~d — > n~pA reaction with 
stopped fT-mesons Here it is the shape of the Ap effective mass spectrum near 
the kinematic limit which is sensitive to the parameters. In impulse approximation 
the amplitude for this process is proportional to that for K~n — > 7r~A and, if the 
Fermi motion in the target deuteron is neglected, the reaction is purely s-wave with 
no spin-flip. The final Ap system is therefore in the same spin-triplet state as the np 
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pair in the deuteron and the values at and r t so determined, which are also shown 
in Fig. m, are consistent with those obtained from the scattering experiments [Jj, § . 

Phenomenological investigations of the Ap interaction by the Jiilich || and Ni- 
jmegen || groups yield low energy scattering parameters in agreement with the 
results of Fig. p], though it is impossible to quantify the systematic uncertainties 
inherent in such models. It should, however, be noted that their spin-singlet poten- 
tial is more attractive than the triplet, which is necessary to ensure the correct spin 
assignments of the ground states of and ^He |7[ . 

The final state interaction in the pp — > K + pA reaction at low Ap effective masses 
is also sensitive to the Ap scattering parameters ||. The recent data on this reac- 
tion very close to threshold ||, taken at the COSY-11 magnetic spectrometer [ TO | 
installed at the cooler synchrotron COSY- Jiilich allow us to extract informa- 
tion complementary to that obtained from elastic scattering because it is possible to 
reach lower centre-of-mass momenta. 

In the present paper we aim to constrain the Ap scattering parameters from the 
shapes of the double-differential pp — > K + pA cross section. In section 2 a simplified 
model is outlined to describe the principal interactions between the three outgoing 
particles. Experimental details, including event selection, are discussed in section 3. 
The FSI model was already used in ref. |§ to determine the precise beam energy by 
fitting the integrated total cross sections as a function of beam energy with assumed 
values of the Ap input parameters. To avoid biassing the analysis, in the present 
work we determine the Ap parameters using only the structure of the Dalitz plot 
at each energy and not its normalisation. For this purpose we apply the maximum 
likelihood method to obtain a map of the confidence levels for the Ap scattering 
parameters. This formalism and the definition of the event weights are shown in 
section 4. Even taking the Ap spin-triplet and singlet to be identical, the resulting 
values of the average scattering length and effective range presented in section 5 
are strongly and systematically correlated in the fit, such that it is hard to quote 
error bars. However, the position of the nearby pole in the Ap scattering amplitude, 
corresponding to a virtual bound state of the system, is much more stable, being 
unbound by (7.7l|;o) MeV. Our conclusions are presented in section 6. 
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2 Model for the FSI in the pp — ► pK + A reaction 

As discussed in ref. ||, if the basic production mechanism is of short range then the 
energy dependence of the pp — > pK + A cross section close to threshold is dominated 
by the available three-body phase-space dp(3), modified by final-state interactions. 
In principle one should consider FSI's in all the three two-body subsystems, A-p, 
A-K, and p-K. Since the strong interaction in the first case appears to be more 



than an order of magnitude larger than for the other two (Tj, [13]], we concentrate on 
the dominant factor f FSI (q) in the A-p system. In addition, however, the Coulomb 
repulsion in the proton-kaon subsystem f c {q P K)i where q p K denotes the cm. mo- 
mentum in pK subsystem, is also important at the low energies pertaining in our 
experiment. Keeping only these two interactions our ansatz for the production cross 
section is 

do- ~ fc(q PK ) fFsi(q) dpi?,) . (2) 

The S'-wave assumption implicit here is justified for our experiment since, for 
excess energies e < 7 MeV, the maximum momentum in any two-body subsystem is 
below 100 MeV/c. It should be noted that the bubble chamber data [[[], || cover a 
higher range of A-p momenta. 

Choosing as independent variables S p k and S\k, the squares of the effective 
masses in the p-K and A-K subsystems, integration of Eq. (0) over the angular 
variables leads to a number of events per pixel in the Dalitz plot distribution of the 
form 



cPa 

dS P K d S\ / N - ./..nsi. • 



/ fc(q PK ) f FSI (q)dp(3) . (3) 

J anelcs 



Structure in the Dalitz plot must be associated with the functions Jfsi and f c . 
Now the Coulomb distortion factor defined in Ref. 0] contains no free parameters 



14H - In contrast the fpsi factor depends on the scattering lengths and effective 
ranges in the triplet and singlet states. Since the spin dependence is expected to 
be small [^|, []], as shown by some of the extracted numbers in Fig. [I], and our 
experiment is not sensitive to singlet /triplet differences, we used mean values of the 
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scattering length a and effective range r in the parametrization of Jfsi- in analogy 



with Eq. (1), we take the popular Watson form for the final state interaction [15 



f FSI (q, a, r) = — — ^ 2 ■ ( 1) 

a 2 q 2 + ( — 1 + ^raq 2 



2' 

A typical Dalitz plot distribution calculated from Eqs. (||) and (|4|) with a = —1.6 fm 
and f = 2.3 fm is shown in Fig. 0. 

3 Experiment 

The measurement of the pp — > pK + A reaction was performed at the COSY-Julich 
synchrotron, using the internal target facility COSY-11 [pTOfl . Outgoing protons and 
positively charged kaons were identified by means of particle momentum reconstruc- 
tion in the magnetic field combined with time-of-flight measurement. The four- 
momentum, and hence the missing mass (MM), corresponding to the unobserved 
A-hyperon was calculated from energy-momentum conservation. Details of the ex- 



perimental technique are given elsewhere [16 



An example of the missing mass spectrum is shown in Fig. |3|. The function 
G(MM) used to fit the peak corresponding to good pp — > pK + A events, is combined 
with the smooth background B(MM). Only events which deviate by less than ±2(7 
from the central value were accepted and a weight w = G(MM)/[G(MM) + B(MM)] 
was assigned to each of them to describe the probability that the particular event 
resulted from a pK + A final state rather than being a background signal. Our data 
were well described by a Gaussian form for G(MM) with o = 0.5 MeV. The closer 
MM is to the known A mass m\, the larger the weight w for the event. Thus w can 
be interpreted as a penalty factor which is imposed on events where the observables 
are significantly modified by the application of the kinematic fit procedure. 

The experimental Dalitz plot for e = 4.7 MeV, shown in Fig. f|a, demonstrates 
that the whole kinematically allowed region is occupied by data, so that there are 
no forbidden zones in our acceptance. This experimental acceptance, calculated 
via a Monte-Carlo simulation with a pure phase-space generator, is however non- 
uniform with a Colosseum-like pattern, as illustrated in Fig. f|b. Comparison of the 
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experimental Daiitz plot (Fig. |]a) with the model calculation (Fig. Q), folded with 
the detector acceptance (Fig. £§b), allows a determination of the model parameters 
a and f at each excess energy e. Although they influence the statistical confidence, 
the relative counting rates at different e are not used, which is important since 
the variation of cross section with energy has already been employed in fixing the 
absolute beam energy [|9] . 

4 Fitting procedure 

The maximum likelihood method was applied to determine best values of the a and 
f parameters, though it must be stressed that these will be strongly correlated due 
to the form of the fpsi factor of Eq. (4). A set of model Daiitz plots Da^(S P K, Sak), 
as defined in Eq. (^|) and (|J), was generated at fixed excess energy e over a grid 
in (a, f) with a G [—6, 6] fm and f £ [0, 10] fm and a step-size of 0.2 fm in each 
variable. 

Defining for brevity x = S p k and y = S\k, the probability V that an event will 
be detected with some (x, y) value is 



V(x,y) = Da,f{x,y) A{x,y) / J D a ,r(x,y) A(x,y) dx dy , (5) 

Daiitz plot 

where A(x, y) is the acceptance function. 

Applying the experimental event weight Wi, the likelihood of a single event is 

h=[V{x i , Vi )p = V? 
and the corresponding global likelihood function 

N 

C(a,f) = Y[VT. (6) 
i=i 

Since Wi is the probability that the i-th event is a pK + A and not a background 
reaction, if the whole experiment were repeated M times with perfect background 
subtraction, then each i-th event would appear Mwi times. The likelihood cor- 
responding to each (xi,yi) point is given by I* = V? IWl and the global likelihood 

N 

function is: C* = U (Vi) Mw > = {C) M . 
i=i 



Both likelihood functions C(d,f) and C*(d,f) have their extrema at the same 
(do, fo) points and thus it is sufficient to perform the experiment once, searching for 
the maximum of the C(d,f) function. 

In order to amalgamate the likelihood functions from measurements at different 
values of e, the M was chosen in such a way that M w % was equal to the number 
of measured events N. 

As shown by Eadie HI?] , the quantity — 2 In [C(d, f)/C(a , f )] has an asymptotic 



X 2 distribution corresponding to the two degrees of freedom, i.e. two parameters 
(a, f). The equation 

In C(a, f) = In C(a Q , f ) - |x|(2) (7) 

defines two-dimensional contours in the (a, f) plane for the desired confidence level 
CL = l-(3. 

5 Results 

Contours in the (a, f) plane of the global likelihood function calculated from a sample 
of about 2400 events measured at six excess energies from 2.7 to 6.7 MeV are shown 
in Fig. [5[ Those solutions with a positive value of a are to be excluded since they 
would imply the existence of a bound hp system and the hyperdeuteron has never 
been found. The branch with negative a yields a long narrow ridge with a very strong 
correlation between the scattering length and effective range, such that it is only a 
combination of the two parameters which is well determined by the experiment. The 
averaged values from the literature [|TJ~|| lie close to the ridge but generally slightly 
outside the 99%-confidence contour. 

The parameter which is in fact well determined by this experiment is the energy 
of the nearby pole in the Ap scattering amplitude. To see this, rewrite the final- 
state-interaction factor of Eq. (4) in the form 

fFSi(q;a,f) - " r "' 5 



(q 2 + llW + Tl) 



The Ap scattering amplitude has poles at q = i^ n , where 
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It is straightforward from Eq. (H) to transform the likelihood function into the new 
variables and the contours in the resulting £(71, 72) are shown in Fig. ||| together with 
the literature values Note that the functional form of Eq. @ clearly shows 

that the data are sensitive only to the magnitudes of the 7, and that the results must 
be symmetric under the interchange 71 ^72- We therefore establish the convention 
that 1 72 1 > |7i|. If 71 is positive then that would correspond to a bound state of the 
Ap system, whereas if it is negative then it is an antibound or virtual state of the 
kind with which one is familiar from the low energy proton-proton singlet S-waye. 
The two branches seen in the maximum likelihood contours of Fig. ||] correspond to 
the mere reversal of the sign of 71. Since it is highly unlikely that there would be 
another singularity of the A-p amplitude very close to zero energy, we can assume 
that I72I 3> I71I, in which case we deduce from our fit that 71 = (— 0.451qj 5 ) fm^ 1 . 



It has been argued jl8] that the single pole limit of letting 72 — > 00 in fact provides 
a better representation of fpsi than the scattering length - effective range form of 
Eq. (||). Our value of 71 corresponds to the average A-p system being unbound by 
an amount (7.7lg;o) MeV. 

In order to obtain values of a and f separately, we must use extra experimental 
information such as for example A-p elastic scattering cross section data. It is seen 
from Eq. (1) that the normalisation of this at zero energy is proportional to the 
square of the scattering length. The low energy data of Ref. |l|, [| are shown in 
Fig. [7] together with a fit to the data on the basis of Eq. (1), where the parameters 
a and f are constrained to lie on the maximum likelihood ridge of Fig. |6|. This is 
achieved with a = —2.0 fm and f = 1.0 fm, corresponding to 71 = —0.41 fm -1 and 
72 = 2.4 fm -1 . The fit therefore confirms that 72 ^> 71 and that the Ap scattering 
data do not realistically allow for a separation between singlet and triplet parameters. 
It is worth noting that our production data are sensitive to much lower values of 
e, as indicated by the arrow, than the scattering data, and it is this region which 
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determines best the value of 71 and hence the position of the pole of the virtual 
bound state. 

6 Conclusions 

Through an analysis of the two-dimensional structure of the Dalitz plot for the 
pp — ► pK + A reaction at fixed energies within a simple final-state-interaction model, 
we have established a strong constraint between the spin-averaged Ap scattering 
length and effective range. The data allow us to fit accurately the position of the 
spin- average virtual bound state 71. Since the data were taken at excess energies 
which are inaccessible to low energy elastic scattering experiments, the results are 
complementary and it is appropriate to make a combined fit of the whole data set, 
leading to new values of (a, f). 

The total cross section data of Ref. || have recently been analysed to determine 
values of (a, f) [[19] . Their argument is, however, somewhat cyclic since a final-state 
interaction with fixed scattering length and effective range was already used in the 
experimental analysis to fix the beam energy ||. In principle therefore these values 
should then be found by the fitting procedure, though there are still of course the 
ambiguities discussed in this paper. We avoid falling into this trap by not using 
the relative normalisation of the event rate as a function of the beam energy within 
our fitting procedure. It is rather the structure of the two-dimensional Dalitz plots 
which fixes our parameter values. 

Though we have implicitly assumed that the Ap system produced from the near- 
threshold pp —> pK + A reaction is the same 3:1 spin-average seen in bubble chamber 
scattering experiments, this is not guaranteed. Just as in the .^"-capture experiment 
[fl, the basic reaction mechanism could favour the production of a particular spin 
combination in the final state. If the present experiment were extended through the 
use of a polarised beam and target, then it would be possible to repeat the current 
analysis separately in the singlet and triplet final states to separate these important 
quantities. 
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Figure 1: A-p scattering parameters for the singlet and triplet states marked with 
open and closed symbols, respectively. Values obtained from experimental A-p elastic 
scattering data of Refs. [1] (stars) and [2] (triangles) are shown as well as the triplet 
values obtained from a K~ capture experiment [4] (diamond). It is only in this 
latter case that an attempt was made to quote errors which are, however, strongly 
correlated. Points deduced from the phenomenological potential models of Refs. [5] 
with solution A (circles) and solution B (circles with additional outer circle) and [6] 
(squares) are also shown. 
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Figure 2: Monte-Carlo Dalitz plot for the pp — » pK + A reaction at an excess energy 
e = 4.7 MeV. The numbers 1, 2, and 3 mark regions of diminishing relative energy 
in the Ap, K + p, and K + A two-body systems, respectively. 
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Figure 3: Experimental missing mass distribution for the pK + subsystem from 
the pp — > pK + X reaction. The fitted peak corresponds to the A-particle from the 
pp —>■ pK + A reaction at an excess energy of e — 4.7 MeV. The grey area is an 
estimate of the background. Only events within the ±2cr band (dashed lines) were 
accepted in the final analysis. 
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Figure 4: a) Experimental Dalitz plot at e — 4.7 MeV containing 776 events; b) The 
acceptance of the COSY-11 apparatus shows a typical Colosseum-like structure. 
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Figure 5: The logarithm of the global likelihood function £(a,f) in the (a, f) plane 
obtained using all COSY- 11 data measured at e between 2.7 and 6.7 MeV. Ex- 
perimental and theoretical values of these parameters are displayed using the same 
symbols as in Fig. |l|. Combining our data with the elastic scattering data of Ref. [1,2] 
leads to the open cross at a = —2.0 fm and f = 1.0 fm. 



16 




CL 70% 
CL 90% 
CL 99% 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 

I y 2 \ (fm -1 ) 

Figure 6: The global likelihood function £(71, 72) transformed into the new variables 
(71, 72). The literature values [1-6] are shown using the same conventions as in Fig. |l[ 
Since the fit function is symmetric in (7!, 72) values, only the curves with 72 > 71 are 
significant, and the data are insensitive to the sign of the j n . The cross at (0.41, 2.4) 
results from combining our data with that of low energy elastic Ap scattering [1,2]. 
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Figure 7: Low energy Ap elastic scattering cross section as measured in bubble 
chambers [1,2] as a function of the cm. energy e. The arrow shows the range of 
energies covered by the present COSY-11 measurement. The combined fit of our 
data with the scattering results lead to 71 = —0.41 fm -1 and 72 = 2.4 fm -1 and this 
is shown here as the dashed line. 
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